Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

Background Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm. Results We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses. Conclusion The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab. Supplementary Information The online version contains supplementary material available at 10.1186/s12863-024-01217-x.

RNAi genes perform different biological activities in plants by controlling microbial (RNA and DNA viruses, viroids, insects, and, more recently, fungal diseases) stresses, expression of protein coding genes and chromatin condensation into heterochromatin [27][28][29][30][31][32][33][34][35][36].Most plants contain four DCLs, with DCL1 causing microRNA synthesis and DCL2, DCL3, and DCL4 essential for 22-, 24-, and 21-nucleotide small interfering RNA (siRNA) biogenesis, respectively [37].The DCL2 and DCL4 have roles in the creation of transacting siRNA, while DCL2 is primarily responsible for the generation of various-sized secondary siRNAs [38].The integrated actions of DCL2, DCL3, and DCL4 play a vital role in disease response and protection [39].The AGO1 is sensitive to virus attacks, confirming the concept that in plants, PTGS represents a defense system against viruses [40].The AGO2 gene is activated in response to the presence of turnip crinkle virus (TCV) and cucumber mosaic virus (CMV) [41].The RDR6 plays a significant role in RNA silencing processes [42] and all RDRs have an impact on plant development and stress responses [43].The exogenously injected VICE12-dsRNA targets the VICE12 gene of apple scab fungus to inhibit its growth and conidial spore formation via the RNAi mechanism [44].The transgenic tomato plant has been developed by applying the RNAi mechanism using exogenously injected CopE/TLR6-dsRNAs, resulting in a greater mortality rate for western flower thrips (WFT) than wild-type tomato plants when fed to WFT [45].
The date palm trees are ancient Asian and African trees that are farmed for their delicious, edible and medicinal value [46].One of the most nutritious fruits in the Middle East and North Africa is the date fruit [47].There are more than 100 million date palm trees in the world [48].It is well-known for its significant nutrients, dietary fiber, natural antioxidants and sources of rich bioactive compounds that are useful in the treatment of neurological diseases and cancer [49].Moreover, date fruits are source of carbohydrates, alkaloids, fatty acids (palmitic, linoleic, lauric, and stearic acid), vitamins, carotenoids, flavonoids, polyphenolic compounds, and tannins, as well as other nutrients such as calcium, magnesium, potassium, and phosphorus [50].Date palm harvest in the world has risen from 4.60 to 8.52 million tons from 1994 to 2018 (FAO.2020).Date palm pollen is high in bioactive substances with phytochemical and nutritional properties, which could boost anti-infertility capabilities [51].Date fruit flesh, peel, and pits have anti-mutagenic, hepatoprotective, anti-inflammatory, anti-diabetic, anti-bacterial, anti-viral, anti-fungal, anti-tumor, nephroprotective, heart disease protective, and anti-cancer characteristics [26].Worldwide, up to 50% newly planted date palm trees and fruits are lossing due to the pathogenic infections and other environmental stresses [52][53][54].To improve this situation, RNAi genes may play a vital role according to the literature review.However, so far, there has been no information regarding these major RNAi gene families in the economically important date palm.Therefore, this study aimed to gather extensive information on the major RNAi gene families (DCL, AGO, and RDR ) that may contribute to the growth and development of date palm.

The data source and descriptions
To explore PdRNAi proteins (DCL, AGO, and RDR) from the date palm (Phoenix dactylifera) genome, we have considered its genome/proteome sequences from the National Center for Biotechnology Information (NCBI) database [55,56] with GenBank accession number: GCA_009389715.1 (NCBI taxonomy ID: 42,345, BioSample: SAMN05011615, weblink: https://www.ncbi.nlm.nih.gov/genome/?term=Phoenix+dactylifera).This dataset has been utilized in comparative genomics study to identify male (Y) and female (X) chromosomes in date palm [57].It has also been utilized in transcriptomics and genome-wide association studies [58,59].In this study, this genome dataset has been utilized to explore PdRNAi proteins by BLASTP search with the query sequences of A. thaliana RNAi (AtRNAi) proteins.Total 17 of 20 AtRNAi proteins (4 AtDCL, 10 AtAGO, and 3 AtRDR) sequence has been downloaded from the Arabidopsis Information Resource (TAIR) database that has been developed by P. Lamesch et al. (https://www.arabidopsis.org/) [60].

Integrated bioinformatics analyses
The integrated bioinformatics studies include BLASTP search, multiple sequence alignment (MSA), phylogenetic tree modeling, functional domain analysis, exonintron structure of predicted PdRNAi of date palm, subcellular location, GO analysis, TFs analysis, CAREs analysis and miRNA analysis.

Identification of PdRNAi proteins Exploring PdRNAi proteins by BLASTP serach
The date palm genome from the NCBI database has been analyzed using the basic local alignment search algorithm for proteins (BLASTP) [61] to explore PdR-NAi genes.The protein sequences of RNAi genes have been retrieved using query coverage (≥ 50%) and E-values (0.00) [9,62].However, only the top-scorer aligned sequences have been considered as the sequences of PdR-NAi genes.The genomic length, protein ID, CDS length, and encoded protein length of PdRNAi genes have been retrieved from the NCBI database.The molecular weight, pI, and GRAVY of PdRNAi protein sequences have been calculated using ExPASy [63].

Phylogenetic tree construction among PdRNAi and AtRNAi proteins to fix the names of PdRNAi proteins
The MSAs of the expected PdRNAi protein sequences have been constructed using the Clustal-W [64] method and the MEGAX [65] package.By performing phylogenetic tree analysis on the PdRNAi genes using the Neighbor-joining technique [66], 1,000 bootstrap repetitions [67] have been utilized to confirm the evolutionary connection.The evolutionary distances have been determined using the equal input technique [68].

Conserved domains and motifs analysis of PdRNAi proteins with respect to AtRNAi proteins
Pfam has been used to search conserved domains in protein sequences [69].The TBtools (a Toolkit for Biologists integrating various biological data-handling tools) has been used to display the Pfam's result [70].The conserved domains have been retrieved from the Pfam database to examine the PdRNAi proteins functional domains.We have investigated the conserved motifs in all of the anticipated PdDCL, PdAGO, and PdRDR proteins using the Multiple Expectation Maximization for Motif Elicitation (MEME) webtool [71].The TBtools have been used to display the motif 's result [70].

PdRNAi genes structures with respect to AtRNAi genes
We have been used the Gene Structure Display Server (GSDS 2.0) to determine the gene structure of the anticipated PdRNAi genes [72].The structure of the targeted genes of date palm has been compared with the A. thaliana gene structure through the exon-intron composition.

Sub-cellular localizations of PdRNAi proteins with respect to AtRNAi proteins
The subcellular location of reported PdRNAi proteins in the cell has been investigated by subcellular location analysis.The web tool WoLF PSORT has been used to predict the subcellular location of the targeted proteins [73].We have been used the TBtools program to display the result [70].

Functional enrichment analysis of PdRNAi proteins
The Gene Ontology (GO) analysis has been carried out using an online database PlantTFDB to confirm the participation of PdRNAi proteins in biological processes (BPs) and molecular functions (MFs) terms [74].The p-values has been calculated using Fisher's exact test with Benjamini-Hochberg adjustments.We have considered a GO term with a p-value < 0.05 as statistically significant.

PdRNAi gene regulatory network analysis
In this study, the associated TF families of the PdRNAi genes in date palm has been analyzed using PlantTFcat, a widely used plant transcription factor database [75].The regulatory network and sub-network has been constructed by integrating TFs and PdRNAi genes and analyzed using Cytoscape 3.9.0[76].From the network, we have identified the hub genes and related major hub TFs of PdRNAi based on the degree of connectivity.Ciselement analysis has been done by Plant CAREs [77] database.We have used the Plant miRNA ENcyclopedia (PmiREN) to download date palm mature_miRNA_ expression and mature_miRNA_sequence [78].After that, we have used Plant sRNA target (psRNATarget) [79] to identify the mature_miRNA_expression and mature_miRNA_sequence IDs that corresponded to our predicted PdRNAi genes.The TBtools program has been used to display this data [70].The detailed pipeline of this study is given in Fig. 1.

Identification of PdRNAi proteins
We have identified four DCL, seven AGO, and three RDR proteins from the date palm genome by top-scorer aligned sequences (Table 1).Therefore, we have considered these 14 PdRNAi proteins for further investigation.
For the convenience of presentation, we have denoted P. dactylifera RNAi protein families as PdDCL, PdAGO, and PdRDR and A. thaliana RNAi protein families as AtDCLs, AtAGOs, and AtRDRs.

Charaterization of PdRNAi proteins Conserved domains and motifs of PdRNAi proteins with respect to AtRNAi proteins
In this study, we have compared the conserved domains of PdDCL, PdAGO, and PdRDR with the AtDCL, AtAGO, and AtRDR proteins.From Fig. 3, we have observed that both PdDCLs and AtDCLs shows almost similar conserved domains, including Helicase_C, Dicer_ dimer, PAZ, RNase III, DND1-DSRM, DEAD, dsrm, ResIII domain.Both PdAGOs and AtAGOs have nearly identical conserved domains, such as Piwi, ArgoN, PAZ, ArgoL1, ArgoL2, ArgoMid, and Gly-rich_Ago1 domains.The Gly-rich_Ago1 domain has been found in PdAGO-1 and AtAGO-1.We have found that the PdRDRs and AtRDRs have almost the same conserved domains, including RNA-dependent RNA polymerase (RdRP) domain (Fig. 3).We have selected 10 significant motifs for PdRNAi and AtRNAi proteins by using MEME-suite analysis.Most of the PdDCL and AtDCL proteins have shown almost similar motif distributions.The annotated conserved motifs 1, 3, 4, 5, and 7 represents Ribonuclease 3 domain, motif 2 represents the Helicase-C domain, motif 6 represents the DEAD domain, motif 8 represents the Dicer dimer domain, and motif 10 represents the PAZ domain, of the predicted PdDCL proteins (Fig. S3A).
The annotated conserved motifs 1, 2, 3, 4, 6, 7, and 8 represent the Piwi domain, motif 5 represent the ArgoL1 domain, motif 10 represent the ArgoL2 domain and motif 9 represent the PAZ domain of predicted PdAGO proteins.The annotated conserved motifs 1-10 represent the RdRP domain of projected PdRDR proteins.The MEME studies of conserved motifs in the PdDCL, PdAGO, and PdRDR protein families have revealed comparable functional diversity in date palms.We have also compared the domains and motifs of PdRNAi proteins with OsRNAi proteins.We have perceived that PdDCLs, PdAGOs, and PdRDRs show similar domain and motif patterns according to OsDCLs, OsAGOs, and OsRDRs, respectively (Fig. S2 and S3B).

PdRNAi genes structures with respect to AtRNAi genes
The exon-intron structures of PdRNAi genes have been analyzed with respect to the AtRNAi genes to investigate the similarities of exon-intron between them.We have observed that the exon-intron structures of PdDCLs, PdAGOs, and PdRDRs are almost identical with the AtD-CLs, AtAGOs, and AtRDRs (Fig. 4).According to the gene structure analysis, the PdDCLs have 19-28 exons, PdAGOs have 22-26 exons and PdRDRs have 3-6 exons, which has been nearly identical to AtDCLs AtAGOs and AtRDRs, respectively.The gene PdAGO-7 has three and PdAGO-2,3 gene has four exons.We have also compared the similarities of exon-intron patterns between PdRNAi and OsRNAi genes and found their similarities like AtR-NAi genes (Fig. S4).

Sub-cellular localizations of PdRNAi proteins with respect to AtRNAi proteins
Subcellular localization investigations has been carried out to learn more about the PdRNAi proteins cellular presence.The majority of PdRNAi, as well as AtRNAi proteins, have found in the nucleus, chloroplast, and cytosol, according to sub-cellular localization study (Fig. 5).Some of our PdRNAi proteins have also found in cell vacuoles, cytoskeletal and mitochondria (Fig. 5).In addition, We have also compared the sub-cellular locations between PdRNAi and OsRNAi proteins and noticed their similarities like AtRNAi proteins (Fig. S5).are found in several biological components.In this analysis, nucl-nucleus, cyto-cytosol, chlo-chloroplast, vacu-vacuole, cysk-Cytoskeletal, mito-mitochondria, E.R.-Endoplasmic Reticulum, plas-plastid, pero-peroxisome.

Cis-regulatory factors of PdRNAi genes
The cis-acting regulatory elements (CAREs) study has been performed to investigate how those regulators controlled PdRNAi genes.According to our findings, the regulatory regions of PdRNAi genes have contained lightresponsive (LR), hormone-responsive (HR), and stressresponsive (SR) related motifs (Fig. 7).The LR motifs including 3-AF1 binding site, ACE (ACGT-containing element), AE-box (activation element), Box 4

(involved in light responsiveness), G-box (CACGTG motif ), GATAmotif, GT1-motif (GGTTAA, G G T A A T T, G G T A A A T, GTTAC, TACAGT, and GGTAAA), I-box ( A G A T A T G A T A A A A
), LAMP-element, Sp1 (specificity protein 1), TCT-motif (polypyrimidine initiator), and TCCC-motif have been found in PdRNAi genes (Fig. 7).The significant number of I-box LR motif has found in the expected PdRNAi genes.The cis-acting regulatory elements ABRE (involved in the abscisic acid responsiveness), O2-site (involved in zein metabolism regulation), P-box (gibberellin-responsive element), TCA-element (involved in salicylic acid responsiveness), and TGA-element (auxinresponsive element) has found in the promoter regions of PdRNAi genes, that are involved in the regulation of hormones (HR).Some of our PdRNAi genes are regulated by the cis-acting elements, including TC-rich repeats (involved in defense and stress responsiveness, MBS (MYB binding site involved in drought-inducibility), MYB (involved in drought-inducibility, LTR (involved in low-temperature responsiveness as stress-responsive (SR).We have also compared the cis-elements between PdRNAi and OsRNAi genes and observed their similarities like AtRNAi genes (Fig. S6).
In root, the expression scores of Pda-miR156b, Pda-miR159a, Pda-miR396a, Pda-miR482 are medium.We have also studied the miRNAs and expressions of OsR-NAi and AtRNAi genes to compare with PdRNAi genes.We have found that OsRNAi and AtRNAi have almost similar miRNAs and miRNA expressions according to PdRNAi genes (Figs.S7A, B).

Discussion
The major RNAi genes RDR, DCL, and AGO function as unique short RNAs that regulate a variety of biological processes and pathways by silencing protein-coding genes [80,81].In this study, we have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome guided by AtRNAi proteins Then phylogenetic tree analysis has been performed to investigate the similarities of PdRNAi with AtRNAi and OsRNAi proteins.To validate the predicted PdRNAi proteins computationally, we have studied their characteristics and activities with respect to both AtRNAi and OsRNAi proteins.Here we considered OsRNAi proteins, since both of rice (Oryza sativa) and date palm belong to the same monocot commelinids clade (Fig. 2 and Fig. S1) [19].In order to investigate the roles of PdRNAi proteins in the protein coding gene silencing, stress response and disease resistance, we have analyzed domains, structures, locations, molecular functions and regulatory factors of PdRNAi genes that are discussed below.
The domains and motifs of predicted PdRNAi proteins have displayed similar patterns compared with AtRNAi and OsRNAi proteins (Figs. 3 and 4 and Figs.S2, S3).All PdDCLs and three PdAGO proteins might be functions as protein coding gene silencers, since all PdDCLs and three PdAGOs belong to the same group of AtDCLs and AtAGOs (AtAGO-3, AtAGO-4, AtAGO-6, and AtAGO-9) protein [82][83][84][85][86][87][88].The PdDCL-1 protein may produce 21 and 22nt sRNAs from small stem-loop RNA substrates, whereas PdDCL-2, PdDCL-3, and PdDCL-4 proteins may produce 22nt, 24nt, and 21nt sRNAs from large dsRNA, respectively [89].The PAZ and RNase III domains in Dicer (PdDCLs) synthesize silencing-associated sRNAs from a larger dsRNA [12,90].The RDR (PdRDRs) synthesize the majority of long dsRNAs in plants [91].The PdRDR-6 protein may be required for the production of trans-acting siRNAs (tasiRNAs) because it cleavages the dsRNA.The PdDCL-4 protein may synthesize the dsRNA from the miRNA-cleavage end, resulting in 21-nucleotide siRNA duplexes formation [92][93][94].More specifically, the identified PIWI domain contains structural similarities with ribonuclease-H enzymes, which indicates that it may be responsible for cleaving the target mRNA [95].The predicted DEAD-box and helicase domain may stimulate date palm to produce secondary siRNAs [96,97].The anticipated PAZ and RNase III domains may operate as an indication, creating miRNAs with comparable sizes [98].The predicted PIWI and PAZ may recognize and cleave the target mRNA [99][100][101].Most of the predicted PdRNAi proteins have been significantly enriched in the nucleus, cytosol, and chloroplast (Fig. 5 and Fig. S5), where the first two locations are known to be involved with the occurrence of the RNA silencing mechanism [102,103].The GO enrichment analysis of PdRNAi genes has detected some vital BP and MF terms that are directly involved with the gene silencing mechanisms and those terms are almost identical with AtRNAi and OsRNAi genes (Table 2).Out of 14 predicted PdR-NAi proteins, 8 proteins are involved with the biological process terms as production of siRNA (GO:0030422) Fig. 7 The cis-acting regulatory elements of AtRNAi and predicted PdRNAi genes.The deep to light green color represent the presence degree of that element with the corresponding genes.(A) hormone-responsive cis-elements, (B) light-responsive cis-elements, and (C) stress-responsive cis-elements [104], RNA interference (GO:0016246) [105], and dsRNA fragmentation (GO:0031050) [106].The GO:0030422 term is associated with the 24nt hc-siRNA and 21nt siRNA production pathways [107].Eukaryotic cells use siRNAs (GO:0030422, GO:0035197) [108][109][110][111] to bind and enhance the degradation of specific endogenous mRNAs generated by dsRNA (GO:0031050.GO:0032296) [112][113][114], to decrease protein synthesis at the post-transcriptional stage.The PdRDR proteins are involved with the molecular function term RNA-directed RNA polymerase activity (GO:0003968) [115], which is associated with RNA-mediated gene silencing in plants [116].Three PdDCL proteins are involved with the terms of the molecular function as ribonuclease III activity (GO:0004525) [117][118][119] which specifically cleave dsR-NAs for RNAi mechanisms in plants [90].Two PdRNAi proteins are involved with the molecular function term siRNA binding (GO:0035197), which is essential to the RNAi mechanisms [111].TFs are key regulatory factors that are involved in the activation or deactivation of upstream signaling pathways and bind to specific promoter regions (i.e., cis-acting elements) of their target genes to increase or inhibit the transcriptional rate of those genes [120].Our detected PAZ-Argonaute [121] (similar as OsRNAi) genes, a member of the PTGS family of TFs, may require for RNAi pathway (Fig. 6).The SNF2 is linked with three PdRNAi genes.The SNF2 protein may supply ATP-driven motor components for remodeling systems that may control gene silencing of PdRNAi genes [122].The miRNAs function as critical regulators of eukaryotic gene expression by targeting mRNAs for break-down or translation inhibition [123].PdRNAi and OsRNAi almost similar miRNAs (Figs.S7A, B).Our detected miRNAs have a variety of significant activities, including gene silencing (miR159a, miR482) [124,125].Thus the predicted PdRNAi genes through sRNAs may play a vital role for the protein coding genes silencing.
We found that all PdDCLs and three PdAGOs belong to the same gene groups of AtDCLs and AtAGOs, which have anti-pathogen activity in Arabidopsis and may protect P. dactylifera from pathogens [40,[162][163][164][165][166].The Ribonuclease III (found in PdDCLs) can resist pathogens [167] and the PdRDR proteins may have antiviral activity [168].Date palm antiviral [169] and antifungal [170] defense may rely heavily on RdRP (GO:0003968) activity [171].The dsRNA-specific ribonuclease activity (GO:0032296) [172] may target several date palm RNA viruses [173].Our detected miRNAs in PdRNAi genes may have a variety of significant activities, including resistance to disease (miR167d, miR2118) [174,175].Thus, our predicted PdRNAi genes might be enhanced the date palm growth and development by regulating the respective protein coding genes against disease risks.

Conclusion
In this study, we have identified 4 PdDCL, 7 PdAGO, and 3 PdRDR genes as PdRNAi genes guided by AtRNAi genes through bioinformatics analysis.To gain a better knowledge of PdRNAi corresponding to AtRNAi and OsRNAi gene families, we have investigated their phylogenetic relationship, domain components, genomic structure, subcellular localization, functional annotations, TFs, cis-elements, and miRNAs.The analysis results of this study have indicated that the PdAGO-5, PdAGO-7, and PdAGO-10 genes might be highly involved in growth and development, and PdDCL1, PdAGO1, PdAGO4 and PdRDR5 might be effective against different stresses in date palm.So far, this is the first study to shed light on the major RNAi gene families of date palms.Therefore, the knowledge gathered from this study might be useful resources to enhanced the date palm production by regulating the respective protein coding genes against biotic and abiotic stresses.However, wet-lab experiments are needed to confirm that PdRNAi genes regulate date palm growth and development through the silencing mechanism of the releted protein-coding genes under various stresses.

Fig. 1
Fig.1The outline of this study

Fig. 2
Fig.2The combined phylogenetic tree.In this tree, the PdRNAi (PdDCLs, PdAGOs, and PdRDRs) and the model plant AtRNAi (AtDCLs, AtAGOs, AtRDRs) proteins have been represented by blue and red, respectively.Here DCLs, AGOs, and RDRs families have been represented by orange, purple, and green circles, respectively

Fig. 3
Fig. 3 The conserved domains of AtRNAi and the PdRNAi proteins.Different color represents different conserved domains.Here Piwi indicates Piwi domain, ArgoN indicates N-terminal domain, PAZ indicates PAZ domain, ArgoL2 indicates Argonaute linker 2 domain, ArgoL1 indicates Argonaute linker 1 domain, ArgoMid indicates Mid domain of argonaute, Gly-rich_Ago1 indicates Glycine-rich region of Argonaut, RNase III indicates Ribonuclease III domain, Dicer_dimer indicates Dicer dimerization domain, Helicase_C indicates Helicase conserved C-terminal domain, dsrm indicates Double-stranded RNA binding motif, DEAD indicates DEAD/DEAH box helicase domain, RdRP indicates RNA dependent RNA polymerase, RRM_1 RNA recognition motif 1

Fig. 4
Fig. 4 Structure of AtRNAi and the predicted PdRNAi genes

Fig. 6
Fig.6 The regulatory network among the TFs with PdRNAi, AtRNAi, And OsRNAi genes.The nodes of the network have been colored based on RNAi genes: AtRNAi: blue, OsRNAi: light green, and PdRNAi: orange.The node shape of the network have been based on RNAi genes: DCLs: triangle, AGOs: octagonal, and TFs: round

Table 1
Basic information of predicted DCL, AGO, and RDR protein families of P. dactylifera (Pd) N.B:The protein names, protein ID, chromosomal location, CDS length, protein length (aa) have been collected from NCBI database.The molecular weight, isoelectric point (pI), and grand average of hydropathicity (GRAVY) values have been collected from the ExPASy.Molecular weights (M.W.) have been measured in Daltons (D) and "aa" means amino acid

Table 2
GO enrichment analysis results with the AtRNAi, predicted PdRNAi and OsRNAi proteins